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A novel, protean, topological soliton has recently been shown to emerge in systems of repulsive 
particles in cylindrical geometries, whose statics is described by the number-theoretical objects of 
phyllotaxis. Here we present a minimal and local continuum model that can explain many of the 
features of the phyllotactic soliton, such as locked speed, screw shift, energy transport and, for 
Wigner crystal on a nanotube, charge transport. The treatment is general and should apply to 
other spiraling systems. Unlike e.g. Sine-Gornon-like systems, our solitons can exist between non- 
degenerate structure, imply a power flow through the system, dynamics of the domains it separates; 
we also predict pulses, both static and dynamic. Applications include charge transport in Wigner 
Crystals on nanotubes or A- to B-DNA transitions. 
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INTRODUCTION. 

The topological soliton, the moving domain wall be- 
tween degenerate structures, ubiquitously populates sys- 
tems of discrete symmetry, most notably the Ising 
model, and appears at different scales and in many 
realms of physics pQ: in mechanical or electrical appa- 
rati [U El [3], superconducting Joscphson junctions [4], 
non-perturbative theory of quantum tunneling [5] and 
particle physics (e.g. Yang-Mills monopoles and instan- 
tons jS], sigma model lumps and Skyrmions [7]). 

Recently Nisoli et al. [8. have discovered a novel kind of 
topological soliton in systems of repulsive particles whose 
degenerate statics is dictated by the intricate and fasci- 
nating number-theoretical laws of phyllotaxis. As the 
phenomena they describe are purely geometrical in ori- 
gin, one would expect this kind of "phillotactic" soliton 
to play a role in many different physical systems at dif- 
ferent scales, wherever repulsive particles in cylindrical 
geometries are present [5]. 

Physically, the problem appears to he highly non lo- 
cal, as energy and momentum is not confined inside the 
soliton but flows though it between the rotating domains 
it separates (see below). Yet, in this article we show that 
a minimal, local, continuum model for the phyllotactic 
soliton that subsumes the effect of the kinetic energy 
of the domains into a modification of the potential en- 
ergy in the familiar equivalent newtonian problem [5], 
can correctly predict its speed, the transfer of energy be- 
tween boundaries, the observed shift in screw angle, and 
charge/density variations, as well as classify its rather 
complicated zoology. 

FROM LEVITOV'S TO DYNAMICAL 
PHYLLOTAXIS. 

In the early nineties, while studying vortices in lay- 
ered superconductors, Levitov [5j showed that a system 
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FIG. 1: Top: the dynamical cactus, repulsive particles, here 
dipoles, on the surface of a cylinder, a is the distance between 
consecutive rings, fl and their angular shift. The dashed 
line is the so-called "generative spiral" [10] . Bottom: for 
a dipole-dipole interaction, energy V versus screw angle Q 
(semi-logarithmic plot) for halving values of the parameter 
a/R starting from 0.5 (lowest line), in arbitrary units e. 



of repulsive particles on a long cylinder presents a rich 
degenerate set of ground states; these can be labeled by 
number-theoretical objects, the Farey classes, that arise 
in phyllotaxis, the study of dispositions of leaves on a 
stem and other self-organized structures in botany [TP] , 
While it is still debated whether Levitov's model is the 
best suited to explain botanical phyllotaxis, it certainly 
imported the same mathematical structures to describe 
a variety of physical systems. Already phyllotactic pat- 
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terns had been seen or predicted outside the domain of 
botany [11] . in polypeptide chains [12], cells of Benard 
convection [13] , or vortex lattices in superconductors [9] . 

Dynamical Cactus 

Recently Nisoli et al. [8 have shown numerically and 
experimentally that physical systems whose statics is dic- 
tated by phyllotaxis can access dynamics in ways un- 
known in botany. They have investigated the statics and 
dynamics of a prototypical phyllotactic system, the dy- 
namical cactus of Figure [T] and found that while its stat- 
ics faithfully reproduce phyllotaxis, its dynamics reveal 
new physics beyond pattern formation. In particular, 
the set of linear excitations contains classical rotons and 
maxons; in the nonlinear regime, both dynamical simula- 
tions and experiments revealed a novel highly structured 
family of topological solitons that can change identity 
upon collision, and posses a very rich dynamics [5]. Un- 
like other topological solitons, the phyllotactic one sepa- 
rates domains of different dynamics; these domains can 
store kinetic energy while rotating around the axis of 
the dynamical cactus; energy and momentum flow across 
the soliton. Moreover as those domains are different, so 
the kinks among them have different shape, characteristic 
speed, and behavior under collision. In some simple cases 
the dynamics of the soliton can be modeled heuristically 
via conservation laws, and certain observables — like its 
speed — -can be computed in this way. A more general 
and complete approach would obviously be desirable. 

The dynamical cactus, a simple phyllotactic system 
depicted in Figure [T] consists of repulsive, massive ob- 
jects holonomically constrained on rings rotating around 
a fixed axis: let a — L/N be the distance between con- 
secutive rings for a cylinder of length L and radius R 
containing N objects; Oi is the angular coordinate of the 
i th particle, and u>i = 9i — 9i-\ the angular shift between 
consecutive rings. The total potential energy for unit 
length is then 

where U is any long ranged repulsive interaction that 
makes the sum extensive and well-behaved, such as a 
dipole-dipole or a screened Coulomb. 

It has been shown that stable structures correspond to 
spiraling lattices [5], where each ring is shifted from the 
previous by the same screw angle u>i = il. The energy 
of spiraling structures as a function of fl is reported in 
Figure [T] for different values of a/R: as a/R decreases, 
more and more commensurate spirals become energet- 
ically costly, independently of the repulsive interaction 
used; for every j there is a value of a/R low enough such 
that any commensurate spiral of screw angle 2ni/j with 



i,j relatively prime, becomes a local maximum, as par- 
ticles facing after j rings becomes nearest neighbors in 
the real space. The minima also become more nearly 
degenerate as the density increases, since for angles in- 
commensurate to 7r each particle is embedded in a nearly 
uniform, incommensurately smeared background formed 
by the other particles [5]. 

One-dimensional degenerate systems are entropically 
unstable against domain wall formation [2]. Kinks be- 
tween stable spirals were found numerically and experi- 
mentally [5]. Numerically it was proved that these kinks 
can propagate along the axis of the cylinder. Experi- 
mentally, the dynamical cactus was observed expelling a 
higher energy domain by propagating its kink. 

As discussed in Ref. [5] such axial motion of a kink be- 
tween two domains of different helical angles confronts a 
dilemma: helical phase is unwound from one domain at a 
different rate than it is wound up by the other. Numeri- 
cal simulations show that the moving domain wall solves 
this problem by placing adjacent domains into relative 
rotation (see supplementary materials in Ref. [5] for an 
animation) . 

Detailed dynamical simulations reveal a complex zool- 
ogy of solitons, along with an extraordinarily rich phe- 
nomenology: kinks of different species merge, decay, 
change identity upon collision, and decompose at high 
temperature into a sea of constituent lattice particles. 

CONTINUUM MODEL. 
Lagrangian. 

To describe domain walls between degenerate struc- 
tures of the dynamical cactus we propose a minimal con- 
tinuum model. First in the continuum limit we pro- 
mote the site index i to a continuum variable: ia — > z, 
6i — * 9{z), the screw angle is now u>i — > lu(z) = ad z 9 
(later we will just assume a = 1). A spiral corresponds 
to lu(z) — d z 9(z) ^constant. A kink corresponds, in the 
lu, z diagram, to a transition between two constant val- 
ues. 

We seek a model that is local and returns the right 
energy when the system is in a stable configuration. To 
express the cost in bridging two stable domains we in- 
troduce a generalized rigidity k. We thus propose the 
(linear density of) Lagrangian as 

C=\lP-V{d z 9)-^(dl9) 2 . (2) 

The first term is the kinetic energy (/ = Xa _1 , X being 
the moment of inertia of a ring). The second accounts 
locally for the potential energy, so that the hamiltonian 
from Q returns the right energy for a static (or uniformly 
rotating) stable spiral like those observed experimentally. 
The last term, the lowest order correction allowed both 
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in the field 9 and in its derivatives, conveys the rigidity 
toward spacial variations of the screw angle and ensures 
that a kink always corresponds to a positive energy exci- 
tation. 

The Euler-Lagrangia equation of motion follows di- 
rectly from Eq. ^ as 

Id 2 t 9 = d z V (d z 8) - K d 4 z 6 (3) 

and corresponds to the conservation of the density an- 
gular momentum Id t 9, where the density of current of 
angular momentum is V' (d z 9) — k d z z 9 (V' is the deriva- 
tive of V with respect to its argument and represents the 
local torque). 

In the following we will prove that, whatever the pre- 
cise form of the actual potential V, and the geometric 
parameters of the problem, there exist solitonic solutions 
of Eq. ([3]) between local minima of V. We will show how 
to derive many qualitative and quantitative results by 
reducing the problem to an equivalent one-dimensional 
newtonian equation, as customary for these kind of prob- 
lems [5], yet with a significant twist. 

Equivalent newtonian picture. 

We seek a uniformly translating solution of Eq. ^ 

6 = 6{z- vt) + wt, (4) 

where the constant w accounts for angular rotation in- 
variance, reflects the underlying transverse structure in 
our system (which is lacking in Llein-gordon-like sys- 
tems), and is the source of much of the new physics. By 
defining s = z — vt, we find d 2 9 — v 2 d 2 9/ds 2 , which leads 
to the equivalent newtonian equation [5] 

«^ = -*» (5) 

whose equivalent potential $ is given by 

$ (oj) = ]-v 2 I(w - u)f - V (uj) . (6) 

a) is a constant that is sometimes useful to redefine as 
uj = t/(Iv 2 ) in terms of r, which has the dimension of a 
torque. We will find that in studying static solitons it is 
easier to use t, while for dynamical ones u> is preferable. 

In the equivalent newtonian picture a soliton corre- 
sponds to a trajectory uj(s) of a point particle that starts 
with zero kinetic energy on the top of a maximum of $(cj) 
and, in an infinite amount of "time" s, reaches the neigh- 
boring maximum of equal height, as detailed in Figure [2] 
A similar treatment applies to Klein-Gordon-like solitons 
and tunneling problems [5J Q3] . Yet the presence of an 
extra quadratic term translates in completely different 
physics and accounts for the dynamics of the domains. 
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FIG. 2: (Color online). The equivalent newtonian picture, 
standard for Klein-Gordon-like solitons or for tunneling prob- 
lems [5] : a soliton is a solution of Eq. (J5J) that corresponds to 
a trajectory u(s) (red solid line [light gray]) of a point particle 
that starts and ends in two equally valued maxima (horizontal 
dotted red lines [light gray]) of the equivalent potential energy 
$ (uj) = ^v 2 I(uj — uj) 2 — V(uj) (black solid line). As the kinetic 
energy must be zero on those maxima, it takes the newtonian 
trajectory an infinite amount of time s (i.e. "time" in the new- 
tonian picture, space in the soliton picture) to reach the "top 
of the hill" with asymptotic speed equal to zero [5]. Unlike 
the Kline-Gordon-like case, the presence of the quadratic term 
7}V 2 1 (u> — uj) 2 (dashed line), which accounts for the kinetic en- 
ergy of the domains, slightly shifts the asymptotic value of the 
kinks (red dotted lines) away from stable static structures u>i, 
UJ2 (black dotted lines), the minima of the real potential en- 
ergy V (uj) (dashed line). Instead the soliton connects spirals 
of angle ui\ = wi + 5i, uis = 0J2 + 82, which in general are not 
stable structures of V. The extra quadratic term also allows 
solitons among non-degenerate configurations, and brings the 
speed of the soliton v into the picture. (Everything is in arbi- 
trary units. The soliton is obtained by numerical integration 
of the Eq. g using the energy profile portrayed in the figure.) 

Shift from stability. 

Because of the quadratic term \v 2 1(uj — uj) 2 in the 
expression of <I>, and unlike the Klein-Gordon- like soli- 
ton, our soliton in general does not connect two stable 
structures ui\, LO2, local minima of V; rather it goes from 
U)i = uji + Si to d>2 = u>2 + 82, bridging spirals whose 
screw angle is slightly shifted away from stability, the 
shift increasing with the speed of the kink (Figure [2]) . 

We find then a most interesting result: a moving soli- 
ton cannot connect two stable structures. At most, it can 
bridge a stable spiral uj\ with one shifted away from equi- 
librium uji =0^2 + ^2, if we choose u) = u>i, but only if 
V(u) X ) < V(uj 2 ). 

These shifts away from stability were observed in dy- 
namical simulations (integrating the full equation of mo- 
tion for the discrete system of particles interacting via 
Eq. [lj in the form of precursor waves propagating at the 
speed of sound in front of the soliton, as in Figure [3] 
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They correspond to the physically intuitive fact that, in 
order to move, the kink needs a torque to propel it. 

Clearly, static solitons (v = 0) can bridge stable spi- 
rals if these are degenerate, i.e. , V(uii) = V(oj 2 ). On the 
other hand, if two domains are degenerate, and are con- 
nected by a moving soliton, they both need to be shifted. 

Boundary conditions. 

The conditions of existence of a soliton between the 
asymptotic domains u)\, C02, are 

r &(Qi) = &fa) = o m 

\ = *(u> 2 ). K ' 

Now, from Eq. Q we obtain the expression of the angular 
velocity for a soliton of speed v, 

9(s) = -voj(s) +w, (8) 

which implies that in general, as the soliton travels along 
the axis, one or both of its two domains must be set into 
rotation, and at different angular velocities. As men- 
tioned, this rotation of the domains is indeed observed in 
dynamical simulations [5] as the mechanism with which 
the soliton unwinds and rewinds spirals of different gain- 
ing angles. 

Unlike sine-Gordon-like one-dimensional topological 
solitons, which separate essentially equivalent static do- 
mains and can travel at any subsonic speed |14j . we see 
that phyllotactic domain walls instead separate regions of 
different dynamics: energy and angular momentum flow 
through the topological soliton as it moves, rather than 
being concentrated in it; as mentioned above, this was 
used to show heuristically that its speed v is tightly con- 
trolled by energy-momentum conservation, phase match- 
ing at the interface and boundary conditions [8 . We 
will show how a more precise version of those heuristic 
formulae can be deduced within the framework of our 
continuum model, and in a more general fashion. 

The angular speed of rotation of the domains depends 
on the parameter w, which along with J) is a constant 
to be determined. On the other hand we can show that 
energy conservation constrains w to w. In fact, when 
the two domains are shifted from equilibrium, and thus 
subjected to a torque, energy flows through the bound- 
aries of the system as the domains rotate. By imposing 
no energy accumulation in the kink we fix w in Eq. @ 
in the following way: the power entering the system at 
the asymptotic boundaries can easily be deduced via the 
Noether theorem, to be 

j±oo = - Hm V'(u(a))e(s), (9) 

s — >±oo 

not surprisingly, the torque times the angular speed. By 
requiring j+oo—j-00 an equation for w, ui is found 

w = v(px + u>2 — £>) (10) 



which will come in hand in many practical cases. 



ZOOLOGY. 
Static Kinks. 

Let us explore static kinks first. Equations ^ reduce 

to 

f T = V'(u>i) = V'fa) , . 

\ r(cDi-wa) =V(u>i)-V(Q 2 ). 1 ' 

The first of the two Equations (JTTJ) tell us that r is a 
torque applied at the asymptotic boundaries of the sys- 
tem, while the second shows that no torque is necessary 
if the topological soliton connects degenerate structures, 
as we anticipated. 

Unlike topological solitons of the Klein-Gordon class, 
here static kinks are allowed between non-degenerate do- 
mains, through an applied torque that shifts the two 
structures out of the minima of V. 



Solitons moving between non degenerate domains. 

Unlike the Klcin-Gordon-like Lagrangian, our La- 
grangian is not invariant under Poincare group. Our 
traveling soliton is not just a boost of the static one, 
as it is clear from Figure [2] In particular the shifts of the 
domains from stability increase with the speed v. By ex- 
panding V (w) around its local minima and using Eq. (6) 
one finds easily that small shifts are proportional to the 
square of the speed of the soliton, or 5%, 82 oc v 2 . 

Let us consider a particular case of practical impor- 
tance. Kinks between non-degenerate structures can 
sometimes be found in the experimental settings, because 
of static friction, after having annealed mechanically the 
cactus with free boundaries. If a small perturbation is 
given, enough to overcome that friction, the dynamical 
cactus will expel the higher energy domains by propa- 
gating the kink along the axis [SJ . Figure [3] reports a 
dynamical simulation of this case obtained by integrat- 
ing the full equation of motion for the discrete system of 
particles whose interaction is given in Eq. (JT|). 

Let us show how our continuum model can reproduce 
a soliton between non-degenerate minima and provide 
quantitative predictions. Free boundaries imply j± OQ — 
0: since v =/= 0, V cannot be zero on both boundaries, 
as we saw before, nor can the angular speed of rotation 
of the domains be. The only possible free boundaries 
solution is one that has V = on one side and d t = 
on the other. Hence the soliton must connect, say, 0)2 = 
LJ2 to [b\ — ui + Si, with V(ijJ2) < V(u>i), which fixes 
u) = L02 = tt>2. Equations Q then fix the speed of the 
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soliton as 



2AV~ 
/Aw 2 



(12) 



with AV = V(u>{) — V{u>2), Alo = u>2 — u>i (in this partic- 
ular case UJ2 = 0J2). Eq. |l2"| ) corrects an earlier formula 
found heuristically via energy conservation, which ne- 
glected the shifts [8] . From the angular velocities of rota- 
tion of the domains in the dynamical simulation (Figure [3] 
third panel), a speed of 22.1 s _1 can be computed and 
compared with the value 23.4 s _1 predicted by Eq. ( |12| : 
the small discrepancy is likely due to energy dissipation 
into phonons - a known effect for solitons in discrete sys- 
tems. 



From Eq. (10) we find w — v C02 which together with 
Eq. ^ implies that the region of lower energy rotates 
uniformly, while that of higher energy, which is shifted, 
remains still: since there is no energy flow through the 
boundaries, the kink transforms the potential energy dif- 
ference between the two domains into the kinetic energy 
of rotation - and vice versa, depending on its direction of 
propagation. (Only propagation toward lower potential 
energies were observed experimentally [8].) 

By expanding V around ui, V ~ \Ic\(lo — lui) 2 , where 
ci is the speed of sound in the spiral of angle uj\ , we ob- 
tain an approximate expression for the shift of the screw 
angle in the region of higher energy: 



Aw 



v 2 /c\ 
l-v 2 /c\ 



(13) 



that clearly holds for v -C c\. As the system is prepared 
in a stable, non-shifted configuration of angle uji, a pre- 
cursor in front of the soliton propagates to accommodate 
the shift Si, as it was seen yet not understood in the 
dynamical simulation reported in Figure [3] [8 . Eq. (13) 



applied to the experimentally measured V(ui) employed 
in the simulation, predicts a shift Si — —0.043 rad, 
which fits the results of the simulation well (Figure [3] 
last panel). 

We will not deal here with collisions between solitons. 
Yet, even without knowledge of the precise form of V, a 
few considerations can be made on asymptotic states of 
the collision depicted in Figure [3] the domains beyond 
the approaching pair of kink anti-kink rotate with oppo- 
site angular velocity 6 — vAuj, v given by Eq. (12). After 



collision v changes sign and so must Aw, as the infinitely 
long domains cannot invert angular velocity instanta- 
neously: the emerging asymptotic configuration must be 
that of a pair of different kink anti-kink, connecting the 
old asymptotic domain with a different nearest neighbor, 
hence the collision metamorphosis already discussed in 
Ref. 0. 
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FIG. 3: (Color online). Dynamical simulation for the colli- 
sion/conversion of two solitons, emitted from free boundaries 
of a dynamical cactus. Top panel: screw angle vs. space (in 
number of rings) and time (in seconds); note the elastic wave 
preceding the soliton. Second panel: energy per particle (mj) 
vs. screw angle (rad) for our system (wi = 1.79 rad, UJ2 = 1-43 
rad, u>3 — 2.40 rad); the two kinks - before collision - connects 
a high (inner: u)i) to a low (outer: UJ2) energy domain; after 
collision a low (inner: UJ3) to a high (outer: Third panel: 



angular speed (s 



space at a given time: the speed of 



the soliton can be extracted as v = A6/Alj = 22.1 s _1 (pre- 
dicted 23.4 s~ x ). Fourth Panel: the precursor; plot of the 
screw angle vs. space at different times while the soliton and 
its preceding wave advance. The amplitude <5i of the precur- 
sor is predicted via Eq. [l3]as <5i = tjj\ — Qi = —0.043 rad, in 
excellent agreement with numerical observations. The simu- 
lation uses the experimental density and magnetic interaction 
as in Ref. [8]. (See supplementary materials in Ref. [8] for an 
animation.) 
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FIG. 4: (Color online). Equivalent newtonian diagram (see 
caption of Figure [5] for explanation) for symmetric solitons 
between degenerate spirals. As the quadratic term raises all 
the maxima of — V (uS), solitons with the same speed, but con- 
necting different domain are possible under different applied 
torque. The smallest kink (solid red line [light gray]) and a 
larger kink (solid blue line [dark gray]) are shown. 



Pulses: propagating and frozen-in. 

The equation of motion also predicts pulses, both dy- 
namic and static (frozen-in) . The top panel of Figure [5] 
shows the pulse soliton in the equivalent newtonian pic- 
ture as a trajectory falling from a local maximum of 
say toi, and coming back to it without reaching the neigh- 
boring Cii2 because of a higher potential barrier. 

Clearly pulses posse less inertia than kinks, since only 
the region occupied by the pulse rotates during propaga- 
tion. Theory suggests that these pulses should be able 
to propagate with free boundaries and that no applied 
torque is necessary - although solutions corresponding to 
different boundary conditions, and hence applied torque, 
can be found as well, by choosing w ^ u>i. 

Let us consider pulses in stable structures and thus 
Co = When v > c, the speed of sound in the spiral 
of angle uj\, then oj\ becomes a minimum for $ and we 
obtain a solution that oscillates around lo^. the speed 



Solitons moving between degenerate domains. 

To investigate the degenerate case, typical of dynami- 
cal phyllotaxis of high R/a ratios (Figure fTl, we consider 
for simplicity the potential V of FigureHl degeneracy 
requires now ui = (u>i + 0J2}/2, and if the potential is 
symmetric as in the figure, both boundaries are equally 
shifted in opposite directions and subjected to a torque 
of opposite sign and equal intensity (from Equations ^ ) 



M = ^ 2 (w 2 - wi). 



(14) 



Exactly as in the case of the boosted Klein-Gordon-like 
one, our symmetric soliton becomes shorter as v increases 
(Figure [2), although via a completely different mecha- 
nism. In the Klein- Gordon- like soliton shortening is a 
consequence of relativistic contraction of the Poincarc 
group. Instead, in our case, v raises the difference be- 
tween the maxima of and the minimum among them. 
That means higher kinetic energy in the "valley" of <f> 
for the equivalent trajectory, and thus higher values of 
ckv/ds, hence a shorter soliton. 

Equations ([8j)-( 10 ) tell us that the two domains rotate 
in opposite direction with angular velocity of the same 
intensity 



\6\ = v(lo 2 - u>\). 



(15) 



In practice, the system acts like a mechanical inverter of 
rotation, that transmits a power tQ — Iv 3 Acu 2 /2 along 
the tube. 

For any given velocity, many solitons between degen- 
erate structures are allowed, corresponding to different 
applied torques, as in Figure [4] 
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FIG. 5: (Color online). Equivalent newtonian diagram (see 
caption of Figure [5] for explanation) for pulses. Top: a trav- 
eling pulse among degenerate structures is possible because 
of the non-zero speed v. Bottom: theoretical (red dashed 
[light gray]) and experimental (blue solid [dark grey]) frozen- 
in pulse on the measured energy curve of the dynamical cac- 
tus |S]. 
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of sound is the upper limit for the velocity of the pulse 
soliton. As v decreases three things can happen. 

If V(ivi) < V(oj2), then pulses exist for Vk < v < c, 
Vk given in Eq. ( fl2| ). As v approaches from above the 
critical value Vk, the pulse asymptotically stretches to a 
pair of kink and anti-kink between domains uo\ and 0)2, 
placed at infinite distance from each other, traveling at 
the same speed in the same direction. 

If V(uji) = V(u}2), the situation is the same as above, 
yet the critical value for asymptotic stretch into kink and 
anti-kink is v = 0. 

If V(uji) > V(u>2), we have a pulse soliton for < v < 
Vk, when the velocity is less than the critical value given 
by Eq. (12 1. As v goes to zero the pulse freezes into a 



static one. 

Hence, in absence of torque at the boundaries, our the- 
ory predicts that static pulses can exist in every struc- 
ture, except in the lowest energy one. 

These frozen-in kinks were indeed observed experimen- 
tally (but not understood) in higher energy spirals, per- 
haps arising from a kink-antikink symmetric collision |15j 
at the interplay with friction. Figure [5] bottom, shows 
in solid blue line (dark gray) data from the experimen- 
tal apparatus, along with our frozen-in soliton calculated 
numerically, dotted red line (light gray) for the energy 
$ = V (solid black line) empirically measured in the ex- 
perimental apparatus. 



AXIALLY UNCONSTRAINED CASE. 

We have since now considered only an axially con- 
strained case where the repulsive particles were allowed 
to rotate around an angular coordinate, but not to trans- 
late along the axis. Non-local optimization via a struc- 
tural genetic algorithm has shown that the more gen- 
eral case of axially unconstrained particles on a cylindri- 
cal surface reproduces the same fundamental statics of 
the axially constrained one, i.e. the same spiraling lat- 
tices |15j . Numerical simulations have also shown kinks 
propagating among the spirals [15] . The most significant 
difference with the axially constrained case is a drop in 
density in the region of the kink, to locally relieve the 
mismatch. If these particles were charged, the drop in 
density would endow the soliton with a net charge: the 
phyllotactic soliton could function as a charge carrier. 

In fact, Wigner crystals of electrons on large semicon- 
ducting tubes are candidate environments for the phyl- 
lotactic soliton at nanoscale, as discussed in Ref. [8]. A 
crystal pinned by the corrugation potential and/or impu- 
rities will not slide along the tube under a weak enough 
external field. 

A simple extension of the continuum model can easily 
incorporate variations in the linear density and allow in- 
vestigation of the possibility of charge transport by the 
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FIG. 6: (Color online). Equivalent newtonian diagram for 
both axially contrained and unconstrained soliton. In blue 
dashed (dark gray), the axially constrained soliton among the 
stable structures loi , cu2- The two red solid curves (light gray) 
depict the soliton in the axially unconstrained case-which is 
longer-and its corresponding drop in density A. 



phyllotactic soliton. We rescale the interaction as 



V 



V 



(16) 



where A is the linear density in presence of axial dis- 
placement, whereas A Q = N/L. Let ( = dz/a be the 
relative axial displacement: then, if d z ( <C 1 we have, 
A/A Q = 1/(1 + d z Q). The new Lagrangian reads 



C 



ho 

2 



a2 



2 yi 



v{d z e) 
U + d z C) 2 



2 (W 



(17) 



To gain insight regarding the shape of the soliton we 
restrict ourself to the simpler static case with no applied 



torque. By variation of the lagrangian in Eq. (17) one 
obtains 



d 2 u> 
dz 2 



2V(uj) " 



(i+a,C) 3 



0. 



(18) 



which, along with the normalization condition for A, re- 
turns the density of a static soliton (kink or pulse) as a 
function of co(z) 



A(z) 



V 



V{Q{z)) 



2/3 



(19) 



Here V is the asymptotic value of V at the boundaries. 
The static kink Cu(z) can be found as the the solution of 
the equivalent newtonian Eq. ^ with potential 



<$>{u) = -V^V{oj) 1/ \ 



(20) 



Now, V 1 / 3 and V have the same set of local minima, and 
the same ordering among their values: the extension to 
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axial displacements does not alter any of the conditions 
for existence of kink and pulse solitons, at least for the 
static case. 

Eq. ( 19 1 implies, as expected, a drop in density simi- 



lar to a dark soliton in the region of the kink, and thus, 
for a crystal of charges, a net charge. In particular, the 
higher the potential barrier between the two domains, 
the lower the density at the center of the kink. Also, 
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from Eq. (20 1 we can see that allowing an extra degree 
of freedom makes the kink longer. In practice the kink 
can take advantage of density reduction to relax the po- 
tential barrier between the two domains and can thus be 
longer. All this is detailed in Figure [6] where the soliton 
for both the axially constrained and unconstrained case 
(along its variation in density) is reported for the same 
interaction among particles. Analogous considerations 
apply to static pulses. The case of the axially uncon- 
strained traveling soliton is more complicated and will 
not be treated here. 



CONCLUSION. 

We have introduced a minimal, local, continuum model 
for the phyllotactic soliton and shown that its predictions 
are in excellent agreement with numerical data, that it 
provides a tool to calculate otherwise elusive quantities, 
like charge transport, energy-momentum flow, speed of 
the soliton and angular velocities, and could be used in 
the future to develop its thermodynamics. 

Our continuum model should be applicable to many 
different physical systems. As discussed in Ref. |5], dy- 
namical phyllotaxis is purely geometrical in origin, and 
thus the rich phenomenology of the phyllotactic soliton 
could appear across nearly every field of physics. Indeed 
phyllotactic domain walls have already been seen, but not 
recognized, in simulations [16] of cooled ion beams [17] 
where the system self-organizes into concentric cylindri- 
cal shells. Colloidal particles on a cylindrical substrate 
provide a highly damped version [TS], and polystyrene 
particles in air (as used to investigate [19] the KTHNY 
theory of 2D melting [20]) have reasonably low damping 
and long-range interaction. 

Yet it should be noticed our model has a range of ap- 
plication much wider than pure dynamical phyllotaxis. 
Domain walls in any spiraling system whose energy de- 
pends on the screw angle of the spirals, which manifest 
different stable spiraling structures, and where a gener- 
alized rigidity can be reasonably introduced to describe 
transitions between different spirals, should be described 
by such a formalism. Kinks in spiraling proteins and in 
particular transitions between A and B-DNA might be 
approachable this way. 

The author would like to thank Vincent Crespi and 
Paul Lammert (Penn State University, University Park) 
for useful discussions, and Ryan Kalas and Nicole Jcffcry 
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